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Abstract 

We study the theory of a (global) texture with DBI-like Lagrangian, the higher-dimensional 
generalization of the previously known chiral Born-Infeld theory. This model evades Derrick's 
theorem and enables the existence of solitonic solutions in arbitrary (iV + l)-dimensions. We 
explicitly show the solutions in spherically-symmetric ansatz. These are examples of extended 
topological solitons. We then investigate the coupling of this theory to gravity, and obtain the 
static self-gravitating solitonic p-brane solutions. These non-singular branes can be identified as 
the smooth versions of cosmic p-branes which, in the thin-wall limit, suffers from naked singularities. 
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I. INTRODUCTION 



Topological defects are a subclass of solitons which are static and owe their existence 
to their non-trivial boundary conditions (topology). Examples of defects are domain walls, 
cosmic strings, monopoles, and textures pQ. They are classified based on the homotopy group 
of the corresponding vacuum manifold, In cosmology, topological defects attract 

much attention since they can naturally form due to phase transitions and may have played 
important roles in the early universe J2j 13] . 

Textures are defects that arise when the third homotopy group is non-trivial [U H], 
^(■M) 7^ I. Often called non-singular solitons (due to the phase of the scalar field which 
is everywhere well-defined) they are probably best studied in terms of non-linear sigma 
model [6]. They have been extensively studied in the past due to the possibility of seeding 
galaxy formation [U |5j |8]. Derrick's theorem [9], however, provides an obstacle for the 
existence of stable textures. In (3 + 1) dimensions textures are unstable to collapse. The 
reason is that under spatial rescaling the only minimum static energy configuration is the 
trivial vacuum and thus nothing prevents them from shrinking. One way to evade this no-go 
theorem is by having higher order kinetic term(s) 1 , as has been shown by Skyrme [15], who, 
by adding a kinetic term 4th order in the derivative (called the Skyrme term), obtained static 
solitonic solutions dubbed as Skyrmions. This is possible since the term scales differently 
than the ordinary kinetic term and as a result the defects have a natural scale. Skyrme 
interpreted these Skyrmions as baryons with baryon number coming from the topological 
winding number. 

As the number of spatial dimensions increases Derrick's theorem gives more stringent 
constraint for the defects' stability. In the context of textures the Skyrme term cannot 
support solitonic solutions in higher-than (3 + 1) dimensions. In order to have stable extended 
textures we have to have generalized higher-order (> 4) Skyrme kinetic terms. While it is 
acceptable as an effective field theory the addition of arbitrary non-canonical terms does not 
look aesthetic as they unavoidably introduce new scale(s) in the theory which results in new 
parameters that must be set by hand. A theory with many adjustably-unrelated parameters 
will lose its attractiveness. Moreover, there is no motivation from fundamental physics of why 

1 Recently there has also been several studies on topological defects with non-canonical kinetic term [10HT4] . 
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they should exist. They are simply put by hand to fulfill the purpose of having stable defects. 
On the other hand, the Skyrme term itself is ad hoc. There is no underlying fundamental 
principle dictating such a term should exist, other than to have stable solutions. This led 
people to find generalizations of the Skyrme term that can generate any possible higher order 
terms [T6HT9] . One attractive proposal is the so-called chiral Born-Infeld theory [20J. In this 
theory, the texture appears in the Born-Infeld-type Lagrangian [21] and the higher-order 
terms are generated automatically via Taylor expansion for small mass-scale parameter. 

In recent years field theories with DBI (Dirac-Born-Infeld) kinetic terms have been inten- 
sively studied. They gain interest due to the fact that the appearance of DBI kinetic form 
can be motivated from fundamental theory (i.e., the D-brane Action). They appear in the 
string-theory-inspired inflationary model [22J where the inflaton has DBI kinetic term. In 
topological defects there have been proposals of DBI global strings [23J or (gauge) cosmic 
strings [21]. Gravitational field of DBI global monopoles has also been investigated in [25J. 
All of the discussions focus on the existence of defects in (3 + 1) dimensions. On the other 
hand, global fc-defects in arbitrary (N + 1) dimensions have been shown to exist [2S]- This 
gives an evidence of the existence of stable-static defects with non-canonical kinetic terms 
in higher dimensions. It will particularly be interesting to obtain these extended defects in 
DBI field theories. 

From a completely different point of view, in the last decade theories with extra dimen- 
sions have extensively been studied in the context of braneworld. The so-called braneworld 
scenario suggests that our (3 + 1) dimensional universe is a 3-brane living in a higher- 
dimensional bulk spacetime [27H29] . Vacuum uncharged p-branes solutions in arbitrary D 
spacetime have been found by Gregroy [30]. These solutions are boost-symmetric along the 
branes and asymptotically-flat far from the cores. They, however, suffer from naked singu- 
larities. Since phenomenologically the branes can be modeled by topological defects [3X143^] , 
Gregory argued that an appropriate choice of defects core might smooth out the singular- 
ity and suggested Skyrme model as an example. Indeed, it has recently been shown that 
Einstein-Skyrme model in 7d has self-gravitating solitonic 3-branes solutions in its spec- 
trum [35] . The solutions are regular at the core and and the energy- momentum tensor falls 
fast enough outside so that asymptotically approaches the flat vacuum p-brane metric. This 
is an example of codimension-3 solitonic defects in higher- dimensional theories. To the best 
of our knowledge, so far there has not been any discussion of codimension-higher-than-3 
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solitonic defects that can cure the singularities of the thin-wall p-branes. This is simply 
because there are no solitonic defects (with ordinary kinetic term) in codimension higher 
than three. We might think of using global defects, as in [31], but since their static energy 
is divergent in flat space their coupling to gravity create deficit angle in the metric and thus 
cannot be asymptotically flat [36J. 

In this paper we attempt to answer the following two questions. First, can solitonic defects 
with codimension n > 3, which we dub extended-solitons, exist? Second, can we smooth out 
the singular uncharged p-branes in arbitrary dimensions? Using the generalization of chiral 
Born-Infeld theory [20J we found that the answer for both questions is positive. 



II. REVIEW OF CHIRAL BORN-INFELD THEORY 

Here we will review the theory of chiral Born-Infeld [2H] . As in the Skyrme model [T5] the 
chiral Born-Infeld theory can be conveniently formulated in terms of a chiral field U(t, x), a 
unitary 2x2 scalar matrix transforming under SU(2). The Lagrangian takes the following 
form 2 

C = -P 2 (l - ^l-^Tr(V>)) , (1) 

where L M = Wd^U is the left-chiral current and (3 is a mass dimensional scale parameter 
of the model. The non-trivial topology manifests itself in the existence of the topological 
charge [15] 

B= 2^J d3 ^ kTr ( • (2) 
At (3 — > oo the theory reduces to 

c~--Tt{l,l»), (3) 

the ordinary non-linear sigma model. The advantage of this model is that it does not even 
require the Skyrme term to stabilize the defect. The DBI-form is sufficient to render the 
texture defects stable. 

Considering the spherically-symmetric ansatz 

U = e iFir)f - f , (4) 

2 Here we ignore the / 2 prefactor as it is only meaningful in the context of effective meson field theory. 
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where F(r), the so-called chiral angle, is a function of radial only and r are the Pauli 
matrices, the static energy functional is given by 

POO 

E[F] =8vr/3 2 / (1-R)r 2 dr, (5) 
Jo 

with 



1 / F' 2 sm 2 F\ 

and "primes" denotes derivative with respect to r. The variational principle yields the 
equation of motion 

(7) 



r 2 - , 

R ) R 

The DBI form of the kinetic term ensures that the defects are stable against spatial rescaling 
and using numerical technique the author of [20] explicitly shows the solitonic solutions for 
topological charge B = 1. 



III. EXTENDED CHIRAL BORN-INFELD SOLITONS 



A. The Model 

From the previous discussion we can learn that since the Lagrangian ([I]) can be expanded 
to any order there is hope to stabilize the solitonic solutions in any arbitrary dimension. 
We therefore wish to obtain extended solitons in (N + l)-dimensions. For our purpose, it is 
more convenient to re-write the chiral DBI Lagrangian in terms of non-linear sigma model, 

C-dbi = ftl; 

where i = 1, 2, N + 1. Notice that we express the mass scale parameter as /3n, since it is 
a dimensionful parameter with dimension /3% ~ [M]^ iV+1 ). It is more convenient to make a 
change of variables to the dimensionless units, 

<P^M^<f), P^M^p N , (9) 

where M is the mass scale of the theory. 
For large (3, the Lagrangian reduces to 

C ~ -d^d^ 1 , (10) 
4 



1 + -T^d^d^ 1 - 1 

Pn 



(8) 



the ordinary non-linear sigma model. 

The non-linear constraint = 1 spontaneously breaks the symmetry of the theory, 
SO(N+ 1) — > SO(N). From homotopy theory we know that it defines the vacuum manifold 
M = SO(N + l)/SO(N) = S N , whose iV-th homotopy group is non-trivial, n N (S N ) = Z. 
This topological invariant number manifests itself in the topological current 

.f^n^ = _J_ e fJ.aia 2 ...a N - 1 fo°-2 Q foa 3 Q ^a N (11) 



B. Scaling Argument 

Having examined the non-trivial topology of the vacuum manifold, we are curious whether 
there exists a static defect in its spectrum. Let us consider the static energy 



E = J \1-^1-E 2 j d N x, (12) 

where, for simplicity, we set (3 = 1, E 2 = <9;0 a <9 l a is the gradient energy density 3 . Under 
the scale transformations, the energy rescales as 



E\ = J \- N (1 - y/1 - X^E 2 j d N x. (13) 

Stability implies [9] 

JE 

= 0, 



dE 
~d\ 
d 2 E 



A=l 

> 0, (14) 



A=l 



which translates into 



E 2 



N(l -y/l-E 2 



2 + (1 2N)E 2 + N ^ + _ ^Y^E~ 2 ) > 0, (15) 



(1-E 2 ) 3 / 2 JT^E 2 
respectively. It is immediately clear that there is no stable defect for 4 iV < 3. On the other 
hand, with some algebra one can prove that for each N > 3 there always exists a non-zero 



positive value of E 2 that satisfies condition (15); i.e., the possibility of having static and 



stable defects are not ruled out in any higher dimensions. 



3 < Ei < 1 to give a finite and real energy. 

4 The N = 2 case remains conformally-invariant, as in the ordinary non-linear sigma model case, which we 
do not investigate here. 
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C. Hedgehog Ansatz 



Armed with the existence of stable defects we consider a spherically-symmetric (hedgehog) 
ansatz 5 , i.e., 



$j = (cos#i, sin^x cos 6* 2 , sin 9i sin 6 2 cos #3, . . . , sin^ sin6> 2 sin^ ■ • • sin^), (16) 



with 



6>i = a(r), 

Oj = 0,, j=2,...,A-l, 



(17) 



where a(r) is the chiral angle and 9j are the angular coordinates of an A-sphere. The 
Lagrangian becomes 

£ = /3 2 [v/T^A-l] , (18) 



where 



(19) 



with "primes" denotes derivative with respect to r. This gives the (A + l)-dimensional 
Action 

S = j d N+l X c, 

= tt N -! J dt dr r N - 1 C, (20) 

where Qn-i is the hypersurface area of an (A — l)-sphere. The Least Action Principle yields 
the equation for the chiral angle a(r), 

(A -1) sin 2a 



a'r N ~ l 



0. 



(21) 



.VI -K\ 2r 3 - N VT^K 

It can readily be checked that at weak-coupling limit, j3 — >■ oo, the equation reduces to the 
ordinary (static) texture equation 



„ 2a' sin a 
a" + — w 0. 



(22) 



In this paper we only consider B = 1 topological charge. 
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D. Solitonic DBI Textures in Arbitrary N dimensions 



To obtain unique solutions, eq.(21) should be supplemented with appropriate boundary 



conditions. We impose the following conditions 

a(0) = 7r, 
a(oo) -> 0. 

They are chosen to ensure topological charge B = 1 and fmiteness of the energy. 



(23) 



Eq.(21) cannot be solved in a closed form. We need to employ numerical method to 



obtain the solutions. In this paper we use the shooting method to solve it. To have solitonic 
solutions, we require regularity around the cores. Thus the expansion about r = is expected 
to take the following form 

a(r) = 7T + «i r + 0(r 3 ), (24) 

where the undetermined coefficient ai becomes the shooting parameter, i.e., the adjustable 
initial gradient of the function such that the asymptotic solutions satisfy the boundary 
condition [23] 

a(r) 

3.0 




10 
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FIG. 1: Chiral DBI soliton (with f3 = 1) in flat N = 3. Here we reproduce the result of [20 



We obtain solutions for several number of spatial dimensions iV and the coupling constant 



(3. They are the extended solitons in (N + l)-dimensional chiral DBI theory. In fig{T] we 
reproduce the chiral DBI soliton obtained in [20J . Fig|2] shows solutions for various spatial 
dimensions (up to N = 6). The profiles are all localized around the cores and quickly relax 
to zero outside. 
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FIG. 2: Extended DBI solitons (with (3 = 1/2) in various dimensions (N = 3,4,5, and 6). The 
value of (3 is chosen to emphasize the non-linearity of the theory, where < (3 < 1. 



From figj2] we can see that the higher-dimensional defect profiles tend to fall-off faster 
asymptotically than their lower-dimensional counterparts. This fall-off behavior can easily 



be understood by studying the asymptotic expansion of eq.(21). At large radius, the field 
can be approximated by 



a(r) « ctoo + 5a(r), 
~ 5a(r), 



(25) 



where \5a\ 1 and on the second line we set = due to the boundary condition (23) 
This yields the asymptotic equation 



6a" + (N-l)— «0, 
r 



(26) 



with (non-trivial) solution 



a{r) « 5a{r) oc -^ =I , (27) 



where 6 is some numerical constant. As the dimensions increase the inverse power of the 
polynomial increases, and as a result the chiral angle profile falls off faster. 

For any given number of spatial dimensions, the solution is characterized by the coupling 
constant (3. It is responsible for the size of the defects. We can make some rough estimation 
about the defect core thickness. The defect thickness is the size that minimizes its static 
energy configuration. In domain walls, for example, the size can be estimated by the balance 
of the gradient and potential energy. In our theory, the thickness of defect core is determined 
by the balance of linear gradient energy term and its non-linear counterparts. The energy 
is roughly given by 



£~^ 2 U/l-^-l), (28) 



and so the size can roughly be estimated as 

/T 1 . (29) 



^N{N-2) 

In fig|3] we show solutions for various /3 (in N — 3). For (3 > 1 the theory is weakly- 
coupled and the Lagrangian ^ can be Taylor-expanded to yield the ordinary non-linear 
sigma model with higher-order kinetic correction terms. On the other hand, for < (3 < 1 
we can no longer neglect the higher-order expansion terms. The theory becomes strongly- 
coupled. One can see from fig|3] that the defect size is inversely proportional to /3; as (3 
grows, the defects become thinner. We can understand it from the fact that for large (3 
it effectively reduces to an ordinary texture defect, and texture does not possess a natural 
scale; it is unstable against collapse to trivial vacuum configuration. As (3 — > oo the defect 
is infinitely thin (vacuum). 



IV. DBI TEXTURE BRANES 



Having established the solitonic solutions in higher-dimensional flat case, we would like 
to know whether they can be good models for non-singular p-brane cores. In other words, we 
attempt to prove Gregory's conjecture in [30J that the naked singularity in p-branes vacuum 
solutions can be smoothed-out by an appropriate choice of cores, and the natural candidates 
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FIG. 3: 3 + 1 dimensional DBI solitons for various values of (3. 

would be uncharged topological defects with asymptotically-flat space-time metric. Those 
conditions put very stringent constraints on the choice of defects, and we will show that the 
extended DBI texture satisfies all the requirements. 



A. Cosmic p-Branes in Isotropic Gauge 

Consider D-dimensional 6 p-branes in isotropic gauge, 

ds 2 = B(r)\ v dx"dx v - H{r) 2 {dr 2 + r 2 dn 2 N _ x ) 

where 



with Ti = 1 and 



AT-1 
i=l 

i-1 

Tj>i = sin 2 9i Y[ sin 2 9j, 

i=2 



(30) 



(31) 



(32) 



D = N + p+ 1. 
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and iV is the number of the codimensions of the brane. 
The Ricci tensor components are given by 



Rq 



B" ( B'\ , B'H' , B' 



B" H" H 1 B'H' 

{P + ^ + {N-l) lf + {N-l) irr -{p + l) — 



Ro 



-(N-l 
H" 



IT 



TTI \ z TTI TJl TTI TJl 



(33) 



The vacuum Einstein equation for zeroth component, Rq = 0, can be written as 

[{B p+1 )'H N - 2 r N - 1 ]' = Q. (34) 

By appropriately identifying the constant of integration of equation above, we could try an 
ansatz for B(r), generalizing our isotropic solution for p = N = 3 |35j.: 

B(rY +1 = ( ^— 1 . (35) 



1 + 



'r a \N-2 



with r Q integration constant. This yields 



H(r) 



1 + 



'ro\N-2\ n 



1 - 



'to\N-2 



1 - 



Tn N 2(AT-2)' 



(36) 



which immediately solves R 9 e = 0. The only left parameter to determine, m, can be expressed 



in terms of variables p and iV by plugging (35 ) and (36 ) into the remaining equation, R r r = 0, 
which then, after some algebra, gives 



m 



( p + iKjv-i) 

{N + p-l) 



(37) 



Eqs (35) and (36), supplemented with (|37h , then completely solve the vacuum Einstein 



equations in isotropic form (see also [37]). 



B. The Defect Cores 



The full Action for our model is given by 

5 ' ^ D 



(38) 
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The matter-sector Lagrangian still takes the form of (18), with K now redefined as 



1 fa' 2 (iV-l)sin 2 a; 
= 0* YH 2 + H 2 ^ 2 



The equation for the chiral angle now becomes 



B P+l H N-2 r N-l al 



(N - 1)^+1^-2^-3 gin 2a 



VT^K 

Armed with the energy-momentum tensor 

2 5S 



which yields 



T 



AB 



V9$9 AB 
2 ( -\ 



-qabL + 



1/2 ' 



1 + wdc&dCfa 



(39) 



(40) 



(41) 



rpr 



rpti 



-A 
-c 



-c 



a 



12 



sin 2 a 



H 2 r 2 y/T^K' 
we are ready to solve the coupled Einstein's equations, 



(42) 



(43) 



with k = 16tvGd, along with (40). 



C. Self-gravitating Solitonic DBI p-Branes 

We solve the field equations numerically and show some of the solutions in figs.Q-Q. To 
be solitons the solutions must not be singular at the origin. Therefore we require regularity 
around the cores, 

a(r) ps 7r + a.\ r + 0(r 3 ), 
B(r) ps B + B 2 r 2 + 0(r 4 ), 

H(r) ^ H + H 2 r 2 + 0(r 4 ), (44) 
12 



where a,\, Bq, and Ho are three undetermined constants 7 in terms of which the higher-order 
coefficients of expansion can be expressed, 



B 2 



H 2 



kB Na( - 2f3 2 H^ + 2f3\ 1 - 



2N(N + p 



Naj 



Hq[ -2N P B 2 - BqHq(3 k + B H$P 2 kJ1 



(45) 



2N(N - 1)B 

Notice that regularity requires the expansion of the chiral angle a{r) to consist only of 
odd-power terms and the metric functions, B(r) and H(r), to consist only of even-power 
terms. 

An example of metric profiles is shown in fig. (J6]) . Here we have more input parameters i.e., 
the number of dimensions on the brane (p) and perpendicular to it (N), and free parameters, 
i.e., the gravitational coupling k and chiral coupling (3. In this work we fix p = 3 since we 
are only interested in the most realistic case, the case where we live on a 3-brane. For each 
number of co-dimension iV we have two free parameters, k and /3. 

a(r) 




40 
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FIG. 4: Solitonic DBI 3-brane in (7 + l)-dimensions (N = 4) with f3 = 1/2. 



7 They become the shooting parameters in our numerical method. 
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For each A^-codimension, we investigate the dependence of solutions on (3. The behavior 
resembles its flat space-time counterparts. The coupling constant (3 controls the thickness 
of defects. Here we only show solutions with 0n = 1/2 (strongly-coupled case). While their 
width is determined by (3, their existence depends on k. Variation of k does not change 
the width, but (as in the case of the Skyrme branes [35J) there exists a critical value, K cr u, 
beyond which no solutions exist. The value of K cr u in each co-dimension N depends also 
on (3. For co-dimension N = 4 and /3 = 1/2 we found that K crit ~ 1/10. Different (3 gives 
different n crit . 

3.0 \ 

2.5 ;\ 

2.0 1 \ 
1.5 1 \ 

1.0 ; \ 

0.5 1 \ 



20 40 60 80 100 

FIG. 5: Solitonic DBI 3-brane in (10 + l)-dimensions (N = 7) with = 1/2. 



Unlike the case of Skyrme Branes [35] here we do not find any second branch of solutions 
upon coupling to gravity. We conjecture that this is caused by the nature of DBI form; i.e., 



there is an upper bound value of K (39) above which the Lagrangian becomes imaginary, 
thus non-physical. The square-root form of the Lagrangian suppresses the existence of upper 
branch solutions. 

The metric solutions, fig. ([6]), are indeed the smooth and regular version of (35) and (36). 
To prove this, we employ the same technique as in [35] (see also [381 [39]); that is, we calculate 
the ADM mass [10], 

1 



T, 



ADM 



lim 



2k 1 



dQ 



(N-1)T 



(JV-l)p 



(46) 
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FIG. 6: Profiles of metric functions B(r) (lower) and H(r) (upper) in (10 + l)-dimensions (N = 7) 
with f3 = 1/2. Both asymptote to the flat space. 

which, for our isotropic gauge, gives 

, DM = figv-i) Hm r(iV _ 13 [vB{r)B ^ r) + {N _ 1)H{r)H > {r)] v (47) 

We can identify the single parameter tq that characterizes the p-brane solutions with its 
ADM mass for each given value of codimension. Plotting the vacuum and chiral-DBI branes 
solutions, figsJ^jSJ they asymptotically match at large radius. Note that while the vacuum 
solutions are singular at finite distance from the origin (ro), the numerical ones are regular 
everywhere. They start to deviate at r ~ r . 

V. CONCLUSIONS 

In this paper we have studied extended solitons in arbitrary (N + l)-dimensional non- 
linear sigma model. This is enabled by considering DBI-kinetic term in the Lagrangian. We 
showed that the solutions are stable under spatial rescaling and constructed, using numerical 
techniques, the explicit solutions having spherical symmetry (i.e., hedgehog ansatz). They 
are characterized by the coupling constant (5. For a fixed co-dimension N , the theory is 
weakly-coupled when (3 > 1 and, upon Taylor expansion, can be reduced to the ordinary 
non-linear sigma model which have static (unstable) defects solutions, textures, for (3 — )• oo. 

15 




FIG. 7: Comparison between numerical and thin wall solutions (35) of B(r) for 3-branes (with 
P = 1/2) in (10 + l)-dimensions (N = 7). 
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FIG. 8: Comparison between numerical and thin wall solutions (36) of H(r) for 3-branes (with 
P = 1/2) in (10 + l)-dimensions (N = 7). 

On the other hand, for < p < 1 the theory is now strongly-coupled and the non-linear 
terms becomes significant. We showed some solutions for several numbers of N and values 
of P, and in principle any arbitrary D-dimensional solutions can be found. 

One motivation to study this theory is due to a conjecture put down by Gregory [30] 
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that the naked singularity suffered by static boost-symmetric (uncharged) D-dimensional 
p-branes can be regularized by a judicious choice of core. While topological defects (soli- 
tons) seem to be the natural candidate, the asymptotically-flat condition severely restricts 
the option. It is well-known that neither domain wall [UJ nor (global) string jl2] can grav- 
itationally be static. The metric of a global monopole, on the other hand, is static but 
has a deficit solid angle and thus cannot be asymptotically-flat jl3]. Textures become an 
interesting alternative, but unless some mechanism is introduced to stabilize it under spatial 
rescaling, the status of its static gravitational field is unclear. In [35] it was shown that by 
considering the Skyrmions, as the solitonic textures with non-canonical kinetic terms, we 
numerically proved the conjecture for p = 3 and N = 3. One purpose of this paper is to 
present the completion of the proof for arbitrary N in a particular (string-theory-inspired) 
model. 

The existence of the self-gravitating solitonic p-branes for each dimension depends on 
the value of k, in a sense that there exists K crit beyond which no static solutions found. 
We conjecture that the fate of super-critical configuration is either (i) collapsing to DBI 
black branes, or (ii) non-static (i.e., inflating) DBI branes. At the moment we have yet 
investigated and hope to explore these possibilities in the future. 

The asymptotic flatness of the solutions provide an interesting possibility of realizing 
the Dvali-Gabadadze-Porrati (DGP) [44J braneworld model in higher dimensions. There 
has recently been some studies of DGP-realization in models with bulk topological defects 
fields (see [4"oT - T4"8] ) . Indeed in ref. [49J it was shown that a braneworld model with (hyper- 
)monopole core in arbitrary dimensions can support the existence of metastable massive 
graviton which is essential in the realization of DGP theory. This proposal shed a light on 
the development of the higher-dimensional DGP models. Finding a metastable graviton in 
our theory can be developed along that way. 
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